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I. INTRODUCTION 



Considerable attention has been focused on generalized parton distributions (GPDs) 0, 
0, 0, S S| over the past few years. These functions arise from hadronic matrix elements that 
are non-diagonal in momentum space and contain quark or gluon operators separated by 
a lightlike distance. Traditional parton distributions, which arise in inclusive deep-inelastic 
scattering, involve diagonal hadronic matrix elements of lightlike separated operators. Elas- 
tic form factors, on the other hand, are non-diagonal matrix elements of local current op- 
erators and are probed in exclusive electromagnetic reactions, for example. Thus GPDs 
encompass physics of both inclusive and exclusive processes. At leading twist, these struc- 
ture functions describe the soft physics that enters in a variety of hard-exclusive reactions, 
see the reviews iZe iS, iSi, ilfll] • 

Because lightlike correlations enter into the definitions of GPDs, these distributions have 
a naturally simple decomposition for spacelike processes in terms of light-cone Fock compo- 
nents of the initial and final hadronic states. Employing the light-cone gauge, a represen- 
tation of GPDs in terms of overlaps of light-cone wave functions has been derived jlll . . 
This wave function representation is ideal for physical intuition because of the intrinsic 
link to many-body quantum mechanics. A physical picture emerges when one relates 
GPDs to momentum-dependent distributions of partons in the transverse plane |l3l . 3 3 ■ 
The wave function representation of GPDs has another salient feature: the positivity 
bounds jl^ 12) EM 0) HO, 21 1 are transparent. Additionally there is the possibility 



to elucidate the continuity of GPDs in hght-cone time-ordered perturbation theory [22j. 
Despite these virtues, however, little has been done to understand how the reduction prop- 
erties of GPDs arise within the light-cone wave function representation. For example, lack 
of manifest Lorentz invariance makes the polynomiality of GPD moments seem miraculous. 
Consequently the construction of useful models for GPDs based on light-cone wave functions 
is severely limited. 

In an alternate approach, one formulates the GPDs in terms of underlying double dis- 
tributions (DDs) [16|,|23|. The GPDs are then obtained as projections of the Lorentz in- 
variant DDs. Not surprisingly this formalism explains the polynomiality properties of GPD 
moments. For this reason, the DD formulation naturally has become attractive to model 
builders: DDs are used almost exclusively to model GPDs. In spite of this popularity, calcu- 
lations of model DDs have only recently been pursued 2^ 2^ . A drawback to using models 
based on DDs is that the positivity bounds become obscure. Possible marriage of DDs and 
positivity bounds has been addressed in a particular framework ^]. Overall there is still 
little insight in constructing DDs with dynamical content. Furthermore, as we pointed out, 
connection to li ght -cone wave functions and the intuition of the Fock space expansion is 
difficult at best li^], except in simple model scenarios IE 3^ - 



In this work, we extend the simple scenarios so far pursued to build a DD model for the 
proton. We treat the proton as a bound state of a residual quark and two quarks strongly 
coupled in both the scalar and axial-vector diquark channels. The resulting light-cone wave 
function of the proton has appropriate spin structure: containing correlations where the 
residual quark carries the spin of the proton as well as correlations where the quark and 
diquark are in a relative p-wave. We use this model to derive DDs for the proton. Inclusion 
of the spin structure into double distributions is crucial if one wishes to make contact with 
the spin sum rule for hadrons. While this (in essence two-body) model for the proton is 
crude, model parameters can be tuned so that measured electromagnetic form factors are 



2 



well described at small momentum transfer. The quark distributions would need to be 
matched so that resulting GPDs are suitable for phenomenology. 

This paper has the following organization. First we review our conventions for DDs and 
their relation to GPDs in Sec. IHl We obtain three DD functions for the proton since these 
quantities are directly encountered in our calculations. Next in Sec. IllH we present the 
quark-diquark model under consideration. DDs are calculated in this model in both the 
scalar and axial- vector diquark channels. Relevant identities are gathered in Appendix \K\ 
while the details of the derivation appear in Appendix O Intuition about the model and 
its construction is provided in Appendix^ where the effective light-cone wave function is 
extracted from projecting onto the light-cone. Section IIVI presents phenomenological uses 
for the model. The model is tuned to reproduce the Dirac and Pauli form factors of the 
proton for small momentum transfer. Consequently these simple, model GPDs satisfy known 
constraints, including positivity. A conclusion ends the paper (Sec. |Vj). 



II. DEFINITIONS 

To begin, we set forth our conventions for DDs and their relation to GPDs. Moments 
of DDs appear naturally in the decomposition of twist-two operators' matrix elements that 
are non-diagonal in momentum space; moreover, they provide an elegant explanation of 
the polynomiality property of GPDs. The non-diagonal proton matrix elements of twist- 
two operators can be decomposed in a fully Lorentz covariant fashion in terms of various 
twist-two form factors Ank(t), Bnkif) and Cnkit), namely 

(P',A'|V^(0)7Wz3'^^ ■ ■■zB'^">^(0)|P,A) 



„1 _ _ / A\ 
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A:=0 

-wv(Po^«.(i')i:«(;73i^c7„.wP«...P-.(--) ...(--) . (1) 

where the action of ■ ■ ■ ^ on Lorentz indices produces the symmetric, traceless part of 
the tensor, P is defined to be the average momentum between the initial and final states 
= i(P' + P)^, and A is the momentum transfer A'' = (P' - P)^, with t = A^. T- 
invariance restricts k in the first two sums to be even and odd in the last sum. There are 
three Dirac structures in the above decomposition since in general the twist-two currents are 
not conserved, hence a structure proportional to A'^ is allowed. The decomposition above 
into various form factors is in fact ambiguous. Such difficulties in constructing DDs have 



been addressed in the literature j28|, |31|, |32, |33[, and we find the construction in Eq. 
directly in our calculations. The ambiguity of DDs for spin-i systems has been addressed 
in la? 



The above decomposition can be used to define three double distributions as generating 
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functions for the twist-two form factors 



Bnk{t) 
Cnk(t) 



1 i-i-m 

dp / da(3''-^a^F{(3,a]t), 

-1 J-1+\I3\ 
1 

dp / daP''-''a''K{p,a;t), 

1 ^-l+|/3| 
1 rl-lPl 

dp / daP^'-^a^GiP.a-t). 

-1 "'-1+1/31 



(2) 
(3) 
(4) 



As a consequence of the restriction on k in the sums, the functions F{P, a; t) and K{P, a; t) 
are even in a while G'(/5, a; t) is odd. The a; t) and a; t) DDs are similar in form 
to the functions originally employed in 0]. The difference is due to the third DD, G{P, a; t), 
which incorporates the D-term j3l|, among other things. 

Using the operator product expansion, we can relate the moments in Eq. ^ to matrix el- 
ements of a bilocal operator. By construction, the DD functions appear in the decomposition 
of the light-like separated quark bilinear operator 



1 ni-m 

dp / da e-*/3^-+-A-./2 

1 J~1+\I3\ 



z A A ■ z 

«; t) + o/. ^ K,{P, a- 1) - -^G,{P, a; t) 



2M 



MP), (5) 



where we have appended a flavor subscript q in the relevant places and z^^ is a lightlike 
vector, z"^ = 0. 

Now we define the light-cone correlation function by Fourier transforming with respect 
to the light-cone separation z~ 



1 



M^'^\x,U) = - I dz-e'-^''^-'^' ^' 



(P',A'|^,(-z-/2) 7+^, (z-/2)|P,A). 



(6) 



Above the variable ^, or skewness parameter, is defined by ^ = — A+/(2P+). As is custom- 
arily done, we assume without loss of generality that ^ > 0. The correlation function in 
Eq. © can be written in terms of the two independent GPDs Hq{x,^,t) and Eg{x,C,,t) as 



^+H,{x,^,t) + ^^E,{x,U) 



uxiP). 



(7) 



Unlike the DDs, these GPDs are quantities that enter directly into the amplitude for deeply 
virtual Compton scattering (DVCS), for example. Inserting the DD decomposition Eq. (0) 
into the correlator in Eq. we can express the GPDs as projections of the DDs 



H,{x,^,t) 
E,{x,^,t) 



I i-i-m 

dp / da6{x-p- ^a) [Fg{p, a; t) + ^Gq{p, a; t)] 

■1 J-l+\f3\ 

dp / daS{x- P- ^a) [Kg{p, a; t) + iGq{p, a; t)] 

1 J-1+\I3\ 



(8) 
(9) 
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from which we can view the ^-dependence of GPDs as arising from different shces of Lorentz 
invariant DDs. Due to the symmetry of the DDs with respect to a, the GPDs Hq{x,C,,t) 
and Eq{x,C,,t) are both even functions of the skewness parameter ^. 

The GPD Hg{x,^,t) has an important reduction property. Taking the diagonal hmit of 
the hght-cone correlator Eq. (|7j), we recover the forward parton distributions, 

/l~\x\ 
daFq{x,a;t). (10) 

In DVCS, the relevant current operators produce the charge and flavor structure Ylq^q 
since there are two photons and thus the charge squared weighted GPDs enter in relevant 
physical amplitudes. To discover the relation of GPDs to electromagnetic form factors it is 
advantageous to consider the single photon structure J2q define 

i/(a;,e,t) = 5^e,i/,(x,e,t), (11) 

q 

E{x,i,t) = Y,^Mx,i,t). (12) 

Since Gq{P,a;t) is an odd function of a, we have J^^dp J^-j^^^p^daGq{P,a]t) = and 
consequently the sum rules 



j dxH{x,^,t) = y2eq df3 da Fq{f] , a; t) = Fi{t) , (13) 

J-l q J~l J-1+\P\ 

pl pl ^1-|/3| 

/ dxE{x,^,t) = y2eq dl3 da Kq{(3,a;t) = F2{t), (14) 

J-l q J-l J-l+\l3\ 

which relate the zeroth moments of the GPDs to the Dirac and Pauli form factors, Fi 2{t). 

The GPDs satisfy further constraints arising from the norm on Hilbert space: the pos- 
itivity bounds |l^ 12, El; EM HO, 21 1 . These bounds are particularly important for 



comparing with experiment. Model GPDs which reduce to the experimental quark distri- 
bution in Eq. (fTUj) but violate the positivity bounds should not be considered because one 
knows from the outset that rate estimates predicted by such model GPDs are automatically 
wrong. Violation of the positivity bounds is a signal that the model is inconsistent with 
the underlying field theory. Perhaps surprisingly, such violation occurs frequently in many 
standard hadronic models. Of interest to us are the basic bounds for both the spin-fiip and 
non-flip amplitudes 

e{x-0\M'q^\x,U) \ <\lfq{ ^ ) /. ( ), (15) 

which we use as a stipulation in constructing our model. 

Lastly we need to address the negative range of the DD variable f3. Experimentally and 
diagrammatically (3 is strictly positive and crossing symmetry can be used to relate the 
functions for positive and negative values of (3. To this end, we define two functions for each 
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DD 

F^iP, a; t) = FgiP, a; t) ± F^{P, a; t), (16) 

K^{(3, a; t) = Kg{(3, a; t) ± K^{(3, a; t), (17) 

G^{(3, a; t) = G',(/?, «; t) ± G^(/3, a; t), (18) 

where the antiquark contributions are defined by crossing 

F^{(3,a-t) = -F,{-P,a;t), 
Kq{p,a;t) = -Kg{-P,a;t), 
Gq{f3,a]t) = -Gq{-(3,a;t). 

Thus the plus DDs [F^{P,a;t), Kq{P,a;t), and G^(/5, a;t)] are odd functions of P and 
the minus DDs [F~{P,a;t), K~{P,a]t), and G~{P,a;t)] are even functions. In partonic 
language, the minus DDs correspond to a difference in quark and antiquark DDs (the va- 
lence configuration) while the plus DDs are a sum of quark and anitquark DDs (the total 
configuration). 

By virtue of the above definitions Eqs. (|16 |) -(|18 |) . we can remove the explicit negative-/? 
parts from DDs and GPDs by defining 

H^{x,^,t) = Hq{x,^,t)±Hq{x,U), (19) 
Ef{x,^,t) = Eq{x,^,t)±Eq{x,^,t), (20) 

where the antiquark contributions are defined by crossing, analogous to the DDs above. In 
this form, we can rewrite the reduction relations in a more familiar way 

f^{x) = H^{x,0,0) = f \aF^{x,a;t), (21) 

where the total and valence quark distributions are f^{x) = fq{x) ± fq{x). The zeroth 
moment sum rules are thus 



/ dxH'{x,i,t) = y^eq dp da F-{P,a;t) = Fi{t), (22) 

Jo q Jo J-l+13 

/ dxE-{x,i,t) = y^eq dp daK~{p,a;t) = F2{t), (23) 

Jo q Jo J-l+P 

and depend only on the valence configurations. In our simple valence model for the proton, 
we only address quark configurations. The double distribution variable P as well as the 
momentum fraction x are positive below. The above positivity bounds in Eq. (jisp also hold 
for valence and total distributions and amplitudes constructed from the valence and total 
GPDs. 



III. MODEL DOUBLE DISTRIBUTIONS 

To calculate DDs for the proton, as a first step we use only a simple model consisting 
of two quarks strongly coupled in the scalar and axial- vector diquark channels along with a 
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residual quark. We assume that the residual quark is a free particle of mass m. This model 
can be considered as loosely based on relativistic quark models js^ or on the Nambu-Jona 
Lasinio model of the proton in the static approximation, see, e.g., 0, S^l- We keep full 



Lorentz covariance in order to preserve the polynomiality of the moments Eq. ([Q) . Without 
covariance, once cannot deduce the DDs. 

Currently two-loop calculations, in which the diquark structure is resolved, prevent us 
from obtaining analytic results. Thus for tractability, we treat the scalar and axial-vector 
diquark T-matrices as free particle propagators, writing 

respectively. The proton Bethe-Salpeter vertex for our model is thus 

r(A:, P) = ^ X^'^DiP -k)^ rW(fc, P) + ^ x2 ^""(P -k)® ^i^ik, P), (26) 
where the diquark vertices are direct products of spin and isospin factors 

x(^^ = ^(^75C^)®^(^r2), (27) 



xS = ^('^f^C) ® ^i^r^T,), (28) 



and we do not append propagators for the first and second quarks. For simplicity we choose 
the quark-diquark vertex functions to be point-like, namely 

r(^) = U ® U, (29) 
rS= 757.®r,. (30) 

This choice corresponds to modeling only a subset of the eight possible structures for the pro- 
ton wave function. However in more realistic quark-diquark models, these are the dominant 
contribution, see jH, li^l for a complete discussion. The vertex function also contains 
an overall color anti-symmetrization which we suppress throughout. The conjugate vertex 
is given by T{k, P) = CT{—k, —P)^C^. One could modify the point-like bound-state vertex 
with a form factor as is commonly done in more robust diquark model, but in such models 
the issue of the positivity bounds is hard to address. 

In this model, the axial diquark channel does not contribute to the proton's electromag- 
netic form factors. Thus we can only determine the parameters m and msD by fits to the 
Dirac and Pauli form factors. The parameter m^^) could be tuned by fitting the quark distri- 
butions at some scale, however, we shall pursue a simpler course and set nisD = i^ad = ^d- 
Alternately mA£, could be tuned from neutron form factor data. 

To derive the Fg{P,a;t), Kg{P,a;t) and Gg{P,a;t) DDs, we consider the action of the 
twist-two operator Of^^^^■■■■^^" = '^p^'yi^^iD^^ ■ ■ ■iD^'^^'ipg between non-diagonal proton states. 
To make any progress in calculating DDs, we must use the parton model simplification for 

the gauge covariant derivative: and write D'^ = |((9^ — d^). Working in the 
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P = P- A/2 



P' =P + A/2 



FIG. 1: Impulse approximation to the twist-two matrix elements. Here the twist-two operators 
with momentum insertion are denoted by a cross. The diquark spectator is depicted by a double 
line and the initial- (final-) state proton has momentum P {P'). 



impulse approximation (see Fig. P), we have the contributions from both the scalar and 
axial-vector diquark channels. In the scalar diquark channel, we have 



[fc2 - m2 + is] [{k + A)2 - m2 + ie] [{k ~ P + A/2)2 -ml + ie] 



(31) 



The symmetric, traceless tensor F^'^^ --'^" in Eq. (|!?T|l arises from the momentum space tran- 
scription of the parton model operator O^^^i and is given by 

rMw.../.n _ + A/2)'^i ■■■{k + /\/2Y-\ (32) 

To regulate^ the above expression, we choose to smear the vertex in the following covariant 



^ Alternate schemes using Pauli-Villars subtractions often regulate sucli models and are also attractive from 
the perspective of DDs since Lorentz covariance is maintained. These subtractions, however, generally 
violate the bounds in Eq. H15|l . For example, in the NJL model of the pion with two subtractions 0|, the 
positivity bounds, which were ignored by the authors, are violated for small values of —t. For the case of a 
quark-diquark model regularized via Pauli-Villars subtractions, the violations are more severe and persist 
for all values of ~t due to the mismatch of end-point and crossover behavior. This commonly encountered 
problem is discussed in [2^. For these reasons, we have used the smeared vertex regularization above, 
which we also employed previously for the pion [2^. As shown in one can employ a very similar 
analytic regularization 



r 



r 



(fc + A) 



(33) 



and still maintain the positivity bounds. As with the smeared vertex, this regularization corresponds to 
modification of the active quark's propagators but leaves the spectator diquark propagator untouched. 
There are a number of drawbacks to this regularization when compared to the smeared vertex. For a > 1, 
the model form factors fall off too quickly as a function of —t. Setting a = 1, the proton's electromagnetic 
form factors can only be matched at the 5— 10% level. Additionally the model double distibutions derived 
using the analytic regularization vanish at the boundaries of support, which need not be the case .S^] and 
is not so with the smeared vertex. 
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manner 



A2 A2 
7^^^ . (34) 

Since the NJL model (to which our model bears a resemblance) is non-renormalizable, the 
choice of scheme is incorporated into the dynamics and hence the choice of regularization 
should maintain desired properties. For phenomenological estimates of GPDs, the positivity 
bounds are of crucial importance, and our regularization choice respects these bounds, see 
Appendix ^for details. One can view the regularization A < cxd as mimicking the non-local, 
non-perturbative structure of the twist-two, quark- ant iquark vertex. This choice of smearing 
maintains current conservation but does not respect the Ward-Takahashi identity. ^ Thus the 
normalization of amplitudes in Eq. ()26|) is only approximately preserved. The unregularized 
model is set up so that if the quantities were finite, then the normalization of the Dirac 
form factor -Fi(O) = 1 implies the correct quark content of the proton, = 2, and = 1. 
In the smeared vertex regularization, both of these quantities acquire A dependence. To 
remedy this feature, we keep the scalar diquark contribution normalized so that Fi{0) = 1 
and add a A-dependent prefactor Z{A) to the axial diquark contribution. This factor is then 
adjusted so that A^^^ = 1 and consequently the u-quark distribution is correctly normalized, 
Nu = 2. An additional drawback of the regularization scheme is that the induced quark 
form factors do not become point-like for large momenta. While this is inconsistent with 
asymptotic freedom, the model is meant only for use in the low momentum transfer region 
where the diquark substructure is not resolved. 

The involved numerator structure in Eq. (jHT|) complicates calculation of the symmetric 
traceless part of the tensor, but, this can be calculated directly without recourse to explicitly 
writing out such tensors of rank n. The k^k'^ term of the vector propagator leads to exceed- 
ingly complicated contributions to DDs that are not in the spirit of our simple model. Thus 
we present results including only the g'^'^ term in the main text, leaving the k^k'^ piece to 
Appendix O The tricks employed to obtain these complicated contributions may be useful 
beyond this work. The numerical effects of this diagonal approximation (where the k'^k" 
term is neglected) are expected to be small, especially at low momentum transfer (ssi l39l|. 

In order to compactly write out the DDs, we define the auxiliary functions 

2 1,2 

D''\P,a;t)-' = Pml + ^il-P-a) + jil-P + a)-Pil-P)M^-[il-PY-a^]t/A, (35) 
which is the typical energy denominator in both channels and 

= (X^^^ (CC - ST^t - ^t^T + ^a^t). (36) 

which is the typical regularization prefactor for all DDs in this model. Although m and A 
are not discrete variables, we have employed Kronecker deltas as convenient shorthands for 
fixing the values of a and b. Notice the contraction AabD"'^{P, a; t) is an even function of a. 



^ Since the Ward-Takahashi identities only constrain one generahzed form factor in Eq. ^ (i.e. the electro- 
magnetic coupling) and positivity constrains all generalized form factors non-trivially, we prefer to violate 
the former and preserve the latter. Adding the factor Z{A) to the axial diquark contribution allows us to 
gauge the extent of violation of the Ward-Takahashi identity. Moreover, this factor gives the situation a 
remedy. On the other hand, if the positivity bounds are not respected in the construction of the model, 
it is far from clear how to remedy the violation. 
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Calculation of the DDs for the scalar diquark component of the proton yields 



log D 




m + /5M)2 + [(1 - - 

2M{1- /3){m + /3M) 
4:Ma{m + (3M) 



(37) 

where N represents the overall normalization which is fixed by the charge of the proton 
and D""^ is merely an abbreviation for D°'''{f3, a; t). Despite the behavior of the induced con- 
stituent quark form factors, the asymptotics for Fi{t) and F2{t) that result from Eq. (|H7|) are 
(up to logs) and respectively. The difference can be traced directly to the orbital 
angular momentum content of the scalar diquark light-cone wave function, see Appendix IHl 
In the the axial- vector diquark channel (keeping only the g^^ term), we have the following 
DDs 



/F(")(/?,a;t)\ 



-Z(A)iV(25,,, 











X 


logC°' 1 










Co) 





(2m + /?M)2 - 3m2 - [(1 - fj f 
2MP[2m- {1- P)M] 
-AMa{2ni + (3M) 



(38) 



where Z{A.) is the regularization dependent factor akin to wave function renormalization. 
As commented above, the value of Z{h) is chosen to preserve A^^^ = 1. The full DDs are 
then given by the sum of the scalar and axial- vector pieces, e.g., Fq{(3, a; t) = Fjf\l3^ a; t) + 



IV. PHENOMENOLOGICAL APPLICATIONS 

As commented above, our philosophy is to tune the parameters m and rriD so that pro- 
ton electromagnetic form factor data at low momentum transfer are reproduced. This is 
particularly simple, since the axial diquark does not contribute to these quantities. For 
the electromagnetic form factors of the proton, there is high precision data from Jefferson 
Lab 4^ 3| and a recent global analysis and parametrization of . The Sachs form factors 



are known experimentally to about 2% accuracy in the small momentum transfer regime 
and are given by 

G^(t) = Fi(t) + ^F2(t), (39) 

GMit) = F^{t) + F2{t). (40) 

The form factor data are reasonably fit (to ~ 2 — 5%) by the model for m = 0.445 GeV, 
A = 0.465 GeV, and rriD = 0.720 GeV. For these parameter values the normalization 
A^ = 1.52. In Fi g. El we compare this phenomenological form factor fit to the two pa- 
rameterizations of |45| . As is standard, we plot ratios of electric and magnetic form factors 
to the empirical dipole form factor, namely 
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FIG. 2: Comparison of fits for GE{t) and GM{i) with empirical parameterizations. The ratios 
GE{t)/GD{t) and Guit) / iipGoit) are plotted against -t in GeV^. The cm'ves "Arr:l" and "Arr:2" 
correspond to the parameterizations of Gsit) and G^fit) given in Tables I and II of (4^], respec- 
tively. 




FIG. 3: Quark distributions for the proton model. The u- and d-quark distributions are plotted as 
a function of x. 



where the dipole mass squared is = 0.71 GeV^. 

We can also determine the u- and cZ-quark distributions in our model. Since we do not 
have antiquarks, the plus and minus distributions are identical f^{x) = fg{x). In Fig. |3l 
we plot the u- and ci-quark distributions as a function of x. The distributions are properly 
normalized so that Nu = 2 and Nd = 1. As commented above this normalization requires 
a relative A-dependent factor for the axial-diquark contributions, Z{A = 0.465 GeV) = 
1.09. Again this is required because the regularization scheme we have chosen does not 
preserve the Ward-Takahashi identity. Without the extra factor, the violation is ~ 5-10%. 
Notice that the distributions do not vanish at the end-point x = 0. This is typical of NJL 
type model calculations. The point-like kernel is independent of momentum and hence, 
when one writes down the effective wave equation for the proton wave function, one easily 
deduces that the wave function should be non-zero at both end-points. This is true of the 
unregularized wave equation. The fact that the model distributions vanish at x = 1 is due to 
our choice of regularization. Physically it is thus reasonable to think of the regularization as 



11 





FIG. 4: GPDs for the proton model. The u- and d-quark GPDs H^^{x, ^, t) are plotted as a 
function of x for a few values of ^ at t = —0.1 GeV^. Notice the scale of the d-distributions is half 
that of the u. 



mimicking contributions from higher Fock states. The small-x behavior of the model quark 
distributions, however, compares poorly with experiment. Additionally in Figs. |3] and 13 
we have plotted the H^j^{x, ^, t) and ^, t) GPDs. The figures show the GPDs at fixed 

— t = 0.1 GeV^ for a few values of ^. The value of the u- and d-quark GPDs, Hu^d{x,^,t) 
and Eu^d{x, ^, t), at the crossover x = C, are plotted as a function of ^ in Fig. El for the same 
value of —t = 0.1 GeV^ . The distributions plotted in Figs. EHHl are presumably at some 
low hadronic scale intrinsic to the model. One way to determine the scale of the model is 
to use the evolution equations to evolve empirical parameterizations down to a scale where 
the first few moments of our model distributions agree. This procedure is not unique; many 
models can reproduce the empirical quark distributions at higher scales. Also the use of 
perturbative evolution is questionable at best at low scales. While the evolution kernels 
for GPDs are known at next-to-leading order [i^ 47, il, i^E^I, we caution this approach 
for our simple model. Perturbative evolution cannot generate the non-perturbative small-x 
physics which our model lacks and the small-x physics is crucial for relating to DVCS data. 
The leading-twist DVCS amplitude is convolution of GPDs and a hard scattering kernel that 
emphasizes a region where the final-state wave function is evaluated near the end-point. In 
fact, the imaginary part of the amplitude is directly proportional to an overlap of light-cone 
wave functions, where the final state has a; = 0. In order to be useful for experimental 
comparison, a more pragmatic solution would be to augment the model with a realistic 
parametrization of the valence quark distributions. This could be done by implanting the 
realistic distribution via factorization of the /3-dependence of the DDs j2^. This choice, 
while unrealistic for the DDs, may indeed be less problematic for GPDs since there is ample 



12 








- ^=0.05 


0.(. 
0.3 




- ^=0.05 


h 




■ ■ ^=0.10 






■■^=0.10 


- / / 


\\ 
\\ 


- ^=0.20 - 


^ 0.4 




- ^=0.20 - 




\\ 








/■' ' 


n 




>< 0.3 






/■ / 

-/•"' 
/;/ 


\\ 
\\ 
\\ 




0.2 






/•/ 


\\ 
\\ 






/ ./ - _ 




p 

) 


0.5 




0.1 


' ' •'/ 
) 


1 

0.5 



X X 

FIG. 5: GPDs for the proton model. The u- and d-quark GPDs E^^{x,^,t) are pfotted as a 
function of x for a few values of at t = —0.1 GeV^. Notice the scale of the d-distributions is one 
tenth that of the u. 




FIG. 6: Proton model GPDs at the crossover. The u- and d-quark GPDs, Hu^(i{x,^,t) and 
Eu,d{x,^,t), are plotted as a function of ^ at the crossover x = ^ for t = -0.1 GeV^. For this 
value of t, ^max = 0.33. 



allowance for interplay of x, ^, and t dependence. From the perspective of model building, 
one would be merely satisfying the experimental constraints by matching form factors and 
valence quark distributions; whether or not the resulting GPD models are suitable to make 
contact with data remains to be seen. 
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V. CONCLUSION 



Above we have used a simple quark-diquark model to calculate DDs for the proton. The 
model consists of two quarks strongly coupled in the scalar and axial diquark channels along 
with a residual quark. The simplicity of such a model allows for analytic computation of the 
DD functions which contain appropriate spin structure. The simplicity also allows for us to 
make contact with the light-cone wave functions and the equivalent overlap representation 
of GPDs. This toy model study thus enables a comparison between the physical intuition 
of the light-cone Fock space representation of GPDs and the manifestly Lorentz invariant 
decomposition of DD functions. 

We were careful to choose a regularization scheme that allows for analytic determination 
of the DDs and satisfies both Lorentz invariance and the positivity bounds required of GPDs. 
Our model, although toy-like in nature, is used to match the electromagnetic form factors 
of the proton and can be augmented with realistic valence quark distributions or evolved up 
from its intrinsic scale. 

The scope of our continuing investigation is two fold. One direction is to calculate as 
many DDs in simple scenarios as possible. This will give modelers a better sense of the 
form and behavior of DDs and may assist with empirical parameterizations of GPDs. The 
model used here can be easily extended to calculate double distributions for both the axial 
and tensor GPDs. In another nearly orthogonal direction, one can improve upon the proton 
model used here in order to see how various features of a realistic proton wave function 
manifest themselves in processes like DVCS. Proper treatment of the three-body nature of 
the proton requires two-loop calculations. Suitable analytical and computational tools would 
have to be developed in order to extract double distributions in such a case. The model at 
hand, however, has qualitative similarities to the three-quark structure of the proton and 
should be considered a good starting point. 
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APPENDIX A: BASIC IDENTITIES 

For reference we list identities used in computing DDs above. To calculate numerator 
structures involving the symmetric and traceless tensors, we used generalized Gordon iden- 
tities [s^l of the form 

UyinVuxiP) = ^^A'(P') (2{f,r} + [^F]) u,{P), (Al) 

where F is any Dirac matrix, P^' = + and = (P' - P)^. The usual Gordon 
identity is a special case of Eq. ()A1|) , namely for F = 7^ 

ux'iP')Yux{P) = 7^uy{P') {2P^ + la^^K) uxiP). (A2) 
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We also require the following two cases of the general identity. With 75 = i7°7^7^7^ and 
^0123 ^ _|_]^^ have 

uy{P')ri5 U,{P) = ^Uy{P') (A^75 " B^'^^^a^pP,) Ua(P), (A3) 
UxiP')cJ>^''u,iP) = ^Uy{P') {Y^" - 7^A- + 22£'^'^°^P,7^75) uxiP). (A4) 

APPENDIX B: CALCULATIONS ON THE LIGHT-CONE 

Here we include our conventions for projecting quantities onto the light cone. While 
the development and derivations above rely exclusively on manifest Lorentz invariance, the 
light-cone Fock representation provides transparent physical intuition about our model. The 
light-cone wave functions for the model are admittedly simple and thus provide a useful guide 
to understanding the DDs constructed above. 

For any Lorentz vector a^, we define the light cone coordinates 

a± = ^ (a° ± a^) . (Bl) 

The light-cone spinor ux{k,m) satisfies the Dirac equation — m)u\{k,m) = and is 
explicitly given by (e.g., see |52| ) 

ux{k, m) = ^ (y2k+ + pm + ■ k^) Xx, (B2) 



where /3 = 7°, a = 7^7, and the unit spinors Xx are given by 

X| = i=(l, 0,1,0), 

xt = -L (0,1, 0,-1). 

Using Eq. ()B2|1 . we derive the useful product of spinors of different mass, momentum, and 
spin 

uy{k, m)ux{P, M) = -j^= [5a',a(A;+M + P+m) - A dx'^k^Px " ^^^a)] , (B3) 

yk+P+ 

where the notation ax = + ^Aa^ has been employed with the understanding that spins 
correspond to the signs t= +1 and |= — 1. 

Using the light-cone spinors, one can find the Hg{x, ^, t) and Eq{x, ^, t) GPDs in terms of 
the light-cone, non-diagonal matrix element '^(x, ^, t), namely 

M'/ix, t) = -j^= [(1 - e)H,{x, t) - eE,{x, e, t)] , (b4) 
1 — t, 



M',-\x, t) = -—^^= E,{x, e, t). (B5) 
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These expressions enable an alternate means to derive GPDs for the model considered in 
Sec. mil In this approach, one directly inserts the quark bilocal operator between non- 
diagonal proton states and integrates over the relative light-cone energy k~ at the cost 
of sacrificing manifest Lorentz invariance. This description in terms of light-cone Fock 
components, however, is more intuitive than the DD formulation. Moreover, the light-cone 
energy integration clarifies the positivity properties of our model GPDs. We shall not present 
complete expressions for the GPDs on the light-cone, however, the diagonal overlap region 
X > ^is particularly simple to consider and thus we provide the details for the scalar diquark. 
Analysis in the other region, x < ^, is similar but the expressions are more cumbersome. 

For our model, we can find the lowest Fock component's light-cone wave function by 
projecting the covariant Bethe-Salpeter wave function "^{k, P) onto the null surface = 0, 
namely 

/ ^ i± u 1 f dk'Us,{P ~ k,m) f ^. \_ T 
^Lc(a;,k,gi;si,A) = j^p^ J ^/T^ — (^7 ^TsCj^.^lP - /c, m) 

x^^4^^^^(^'^W(P,M), (B6) 

where x is the fraction of the proton's longitudinal momentum carried by the residual quark 
[x = k~^/P^), and the relative transverse momentum is k^j = h^ — xP^. Above, A labels the 
spin of the proton, whereas the label the spins of the three quarks. We have omitted the 
color and isospin parts of the wave function which are trivial: oc eci,c2,c3(5i,n52,d — 5i,d52,n)53,u 
for the scalar diquark. In the scalar channel, the quark-diquark Bethe-Salpeter wave function 
is 

nk,P) = j^i^^{ D(P-k), (B7) 

and above we have chosen a single Dirac component F*^*^ = —ig^'^^ and the coupling constant g 
has been absorbed into the overall normalization. Our choice of vertex functions corresponds 
to only modeling a subset of the possible three quark light-cone wave functions of the proton 
jH^. Carrying out the projection in Eq. ()B6|) . yields 



Vlc[x, k^eb Si, Aj oc r ITn ; ^2 1 • l^^J 



Xy/1 — X 



The smeared vertex regularization, Eq. (jH^ . generates an effective wave function 
'ijj^{x, k^^i, Si, A) because the spectator diquark's energy pole is unaffected, namely 

^iH., ki,; ... A) oc ^^ sA.,^,(xM^m) + XKM ^ ,3^, 



X\^ — X 

where the function 



M2 - 

X(l—X) X 1- 



^c. = jJ^{ST-6^), (BIO) 
A^ — 

denotes the two terms induced by the regularization scheme. In going from Eq. ()B8|) to 
Eq. ()B9|1 . the effect of regularization is to mimic the contribution from higher Fock states. 
However, there are no true higher Fock components in this model since the interaction kernel 
is instantaneous in light-cone time Focusing on the spin structure, the scalar diquark 
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channel consists of two states: a state where the proton spin is ahgned with the residual 
quark's spin and a state where the quark and diquark are in a relative p-wave. Notice the 
regularization Eq. ()34|) does not alter the spin structure of the wave function. 
The quark distribution function can be obtained from 

fu{x) = J2 [ dk^\ijfc{x,'k^;s„X)\\ (Bll) 

and agrees with the covariant calculation of from its moments in Sec. 11111 Similarly, 

the Dirac and Pauli form factors can be expressed in terms of the effective wave function 
since the light-cone contour integration only encompasses the diquark pole. The expression 
are 

Flit) = Y. I c?k^c?a;V^rc(x,k^ + (l-x)A^;s„A)^^|(x,k^;s„A), (B12) 
F2it) = 11 j c?k^t^a;V'Lc(x,k^ + (l-a;)A^;s„-A)7/>£^(a;,k^;s„A). (B13) 



AAa 

2M 



Also for this reason, we can express the GPDs as simple convolutions in the diagonal overlap 
region x > Defining xi = and X2 = jz^, we have 

x^/jll(^xi,k^-{l-x,)^;s,,Xy (B14) 

with a very similar expression holding for the spin-flip amplitude. In the above form the 
positivity bound is manifest. 

In the axial-vector diquark channel we have the orthogonal amplitude 

- "«3(^, P)ux{P, M). (B15) 



We again have omitted the color and isospin parts of the wave function which are: 
^ £ci,c2,c3[i.Knh,d + h452,u)h,u - 25i,«52,«<53,J- In this channel, the quark-diquark Bethe- 
Salpeter vector wave function is 

^^{k, P) = ( - ^^?(")7.75) D^'^iP - k), (B16) 

and we have chosen a single Dirac component vl^^ = (7 ^"■'71/75 and the coupling constant 
gi"-) = g(^)_ One can carry out the projection in Eq. ()B15|) to find analogous formulas for 
the axial diquark contribution to the quark distributions, form factors and GPDs for x > ^. 
The positivity bounds are again satisfied. 
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APPENDIX C: DERIVATION OF THE DOUBLE DISTRIBUTIONS 



In this Appendix, we detail the calculation of the DDs in the scalar diquark channel and 
comment on the calculation in the axial-vector channel. The crucial steps in the derivation 
hinge upon reducing the numerator by factors present in the denominator or by use of the 
simple identity: 

= 0. (CI) 

In the scalar diquark channel, let us take Eq. (j31|) as our starting point. Denote the 
propagators simply by 2t = (/c — P + — mjj + is, ^ = {k + A)^ — + ie, and (T = 

k"^ — m? + ie. The DDs in the scalar channel can be deduced without reducing factors in the 
numerator. We merely introduce two Feynman parameters {x, y} to cast the denominator 
specifically in the form [a;2t + yQS + (1 — a; — y)^^- One then translates k^ to render the 
integral hyperspherically symmetric via the definition k^ = + (3P^ — (a + l)A'^/2. Here 
(3 = X and a = x + 2y — 1. Carrying out this procedure on Eq. produces 

(C2) 



where the numerator Dirac structure is given by 



+m'^Y + ima'^'A^ + 2m{l + (3P - aA jlf. (C3) 
Using Eq. ()C1|) . as well as the identities in Appendix^ we can cast Eq. ()C2|) in the form 

/■I /-l-zS r 

Aab dp da dH [f - D''\p,a;t)-Y^ ux'iP')M^^uxiP) 



J-l+13 



-i:«^;r3i^''"-'«'^"'--^''"-'(-2) ■••(") ■ (") 



where 



Af^ = |(m + /3M)2-^ + i[(l-/?)2-«2]|7M 

+ia'"' A^{1 -p){m + /5M) - a{m + /3M)A'^. (C5) 

The / integration is then standard and the DDs can be read off simply by using Eqs. (Q) - 
dD) and we arrive at Eq. (jH7|) . 

Calculation in the axial- vector diquark channel is similar, however, the numerator is more 
complicated. In the same units as Eq. (jHTj) . we have the contribution from Fig. ^for the 
axial diquark 



^z(A) / d^k D'^^ik - p) ^ ' ' ^ (ce) 

Q J [k^ - a? + ie][{k + Af - + le] ^ ' 

where the axial- vector propagator, D°'^{k — P), is given in Eq. ()25|) . The terms which result 
from the g""^ structure in the vector propagator can be dealt with straightforwardly after 
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evaluating the contracted gamma matrices. The integrals encountered are then similar to 
those in the scalar channel. 

The terms in the numerator which arise from the second Lorentz structure of the vector 
propagator are more subtle. The Dirac structure of these terms is 

V - f ) [-2m( A; + A/2Y + m V - ima'^'^ A, + + (C7) 

One must keep in mind that the index yU is entangled in the symmetric traceless combination 
pAiMi-- /in given in Eq. ()32|) . The first term in the square brackets of Eq. ()C7|1 can be rewritten 
in the form 

2m{k + A/2Y{k - Pf = 2m{k + A/2)'^[2l + ml], (C8) 

where 2t is the energy denominator of the spectator diquark and when canceled gives rise to 
a contribution oc 6{/3) in G{/3,a;t). The remaining term oc m|, is simple to evaluate. The 
second term in Eq. ()C7|) is similarly easy to evaluate. 

The remaining two terms' contributions to DDs are more involved. One must first reduce 
the quadruple and quintuple products of gamma matrices in these terms. After this proce- 
dure, one is left with terms not yet encountered above. After suitable algebra, these can be 
cast in the form^ 

P-l{l + (3P- aA/2)^f' ■■■{l + (3P- aA/2Y--^\ (Cll) 

where we have shown only the relevant part of the numerator. Appealing to hyperspherical 
symmetry and contracting n lightlike vectors, with the above expression, we have 

-nP.zipp.z-aA. ./2)"- = E kl^^-^^^ ' ^""^ ' "^^^^ ' ^^''^ 

Let us denote the result of the /-integration as f{f3, a; t). We have thus cast terms with the 
numerator Eq. ()C11|) in the form 

d/3 da f{(3, a;t) 
Jo J-i+p 

,^, /?"-fea^piA^---P^"-^-M --) ■■■(--) , (C13) 
dp^^k\{n-k)r \ 2) \ 2) ' ^ ^ 

from which we change the order of integration and then integrate by parts to read off 
contributions to the DDs: -df{P,a;t)/dp, -6{P)f{P,a;t) and 6{p -l + \a\)f{p,a;t). The 



^ The case with an overall prefactor of A • Z proceeds as follows: 

A-l{l + 0P- aA/2)^^ •••(/ + /3P - aA/2)^"-i>, (C9) 
when contracted with lightlike vectors turns into 

-uA-zipP-z-aA- z/2)"-i ^^^g-Yl kU^^^'^ ' ^^"^ ' "Z^)' ■ (dO) 

In this form, we integrate by parts to read off contributions to the DDs, including the surface terms. 
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latter two are surface terms. In particular, the term proportional to S{j3 — l + \a\) contributes 
to the DDs at their boundary of support and need not vanish js^l- 

Adding these contributions to those in Eq. ()38|). we arrive at the full expressions for the 
DDs in the axial diquark channel 



+ 5l3mM+ ^/3[m|j + (5/3 - 3)] + ^(5a^ - - 4/3 - 3) - i[(/3 + a)a^ + {/3 - a)b^]\ 
M^[3m - 2(1 - /3)M] 
-2Ma[3m + (1 - 20)M] J 



+A„i, logD"'' |l + i[/3i(|Q| +/3 - 1) + (1 - f3)S(a+0- 1) + (1 - /3)<5(a - ^ + l)]j ^oj 



2m' 



m^(M^(l - fjy + f [(1 - - a^]) - /3t(M^{l - /3)^ + 2(1 - /3)^]) + M^(l - I3)(2M^P - 2m^ + t)\ 

-D"" I 2m2M2(l - /3)(2 - /3) + mAf(2M2(l - /3)2 + a(a^ - b^) + [(1 - /3)2 - 3a^]t/2) + /3M^(|a^ + 3/3(2 - ;3) - 2) 



-2(1 - /3)M[m(6^ - a^) + a(2M^ - 2m-' M + mt)] 

i[a2 + 6^ - 2m|, + 2M^(1 - 2/3) - (1 - /3)t](-2m^ + 2/3M^ + t[l - 2/3(1 - /3)1)\ 
+ D"" I 2/3M2(6^(a - ;3 + 1) - a^(a + /3 - 1) - (1 - /3)[3(2/3 - l)m^ + 2mjjl) - Mt 



9 Q n 
















+ ct)l 













2A/m - 4M^(3/3 - 1) 




- /3M^ + 13(1 - ri)t\ 



-a{M'' + /3t) 

+ logI3"°[i(c« + /3 - 1) - 5(q - /3 + 1)] I 2aM(2l3M - m) I + log D^^l^dal + /3 - 1) - 5(/3)l I -4Af^/3(l - /3) 

-4mM(l - /3) / V 



A^ 
4m?, 



^(^^ + + (1 - /3)| + (2/3 - fjM^j (A„D"'^ + AjD^'') - |(A„D"*a2 + A^D'^'-b^) + i|!:(A„D"A _ AfcD^'') 



A2 



(1 - 0)M-'(AaD'"^ + A^D'^") 
4(1 - /3)M2(AaD°'^ - AbD^') 

[5(c« + /3 - 1) - 5(o. - /3 + 1)] log +_[i(|(j| + ^ _ 1) _ i(0)]log_-_^ 



1 

2m? 



«(/3)A„b logD" 





4c«A-/(Af - m) 



, '5(/3) ,.2 2,, O 

+ a[A — m ) log — 

2 O 



Am /I 





(C14) 



These cumbersome expressions belie the simplicity of calculating model DDs for the proton. 
Moreover the regularization scheme employed is likely inconsistent between the g^'^ and k'^k" 
terms. A particular signal of trouble is that the ci-quark distribution calculated from the 
above expression becomes negative for small x. For these reasons, we have not presented the 
full results for the axial channel in the main text nor used them in our phenomenological 
analysis. It is interesting to note that the largest contribution from the k'^k'^ piece is at 
small-x and thus observables like the magnetic moment and charge radius are not terribly 
sensitive to the diagonal approximation. 
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